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We construct exact ring soliton-like solutions of the cylindrically symmetric (i.e., radial) Gross- 
Pitaevskii equation with a potential, using the similarity transformation method. Depending on the 
choice of the allowed free functions, the solutions can take the form of stationary dark or bright 
rings whose time dependence is in the phase dynamics only, or oscillating and bouncing solutions, 
related to the second Painleve transcendent. In each case the potential can be chosen to be time- 
independent. 
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I. INTRODUCTION 

Many nonlinear equations of motion allow interesting 
solutions such as bright and dark solitons, which prop- 
agate without dispersion [TJ [2] . More generally, we can 
also consider localised soliton-like stationary solutions, 
which have the appearance of solitons except for some 
aspects of dynamics. They are usually accompanied by 
a spatially inhomogeneous external potential-like term in 
the equation of motion and often called also solitons, or 
solitary waves, or for dark soliton cases, kink-like solu- 
tions [3]. 

The nonlinear Schrodinger equation (NLS), 



(1) 



also known as the Gross-Pitaevskii equation (GPE), is 
of special interest since it describes the relevant physics 
in nonlinear optical systems |4] as well as in degenerate 
quantum gases [51 IS]- For V = and a = —1 it is called 
focusing and allows bright solitons, while for a = +1 
it is called de-focusing and allows dark solitons. These 
solutions tend to work well also when a potential term 
V is added to the nonlinear equation. Equation ([T]) is 
one-dimensional, and while the trivial extension to two 
dimensions by replacing the a;-derivatives by a Lapla- 
cian is possible, the corresponding solutions are plane 
waves without solitonic properties 3'. Typically they 
are not stable: the stationary bright soliton solution sur- 
vives only under certain conditions and if supported by 
an external potential [7j, and the stationary dark soli- 
ton decays into vortices by the snaking instability unless 
the system is in practice a quasi-one-dimensional one . 
This applies especially to any static solution of Eq. (Il|, 
which gives the motivation to seek such multidimensional 
solutions that can be expected to have at least extended 
lifetimes. The static soliton-like solutions are also inter- 
esting more generally, connecting to kinks and domain 
walls [3], although e.g. in superfluid He the multicompo- 
nent structure of the order parameter leads to a rich set 
of complex structures . 

Instead of making such a plane- wave extension into two 
dimensions, one can look at cylindrically symmetric sys- 
tems, which brings forward the concept of ring solitons. 



Since their introduction in a nonlinear optics setting |10j , 
the interest has been in dark ring solitons, which have 
been treated essentially numerically and in the domain 
of a large radius [HI [12] . Although one can reduce the 
dynamics into a one-dimensional radial equation, exact 
analytic solutions (either bright or dark) have not been 
obtained so far. It has been only shown that in the limit 
of infinitesimal amplitude, such solitons are analytically 
described by the radial KdV-equation [TU]. Apart from 
being simply convenient, exact solutions would allow one 
to analyse the stability and decay of ring solitons, includ- 
ing the effect of external potentials, in the same fashion 
as for trivially extended plane-wave solutions in two di- 
mensions. 

We have found that using a similarity transformation 
method, one can actually construct exact analytical ring 
soliton-like solutions of the radial Gross-Pitaevskii equa- 
tion. They are limited to certain forms of external po- 
tential, but even then, they can provide a starting point 
for further studies of the ring solitons in cylindrically 
symmetric potentials. In this paper we first describe the 
similarity transformation in Sec. |ll) and then present the 
solutions in Sec. 



Ill 



Then, in Sec. IV we discuss the 
similarity solutions associated with the second Painleve 
equation. We summarise our work and discuss its impli- 
cations in Sec. (V] 



II. THE SIMILARITY TRANSFORMATION 



The radial GPE in dimensionless form is given by 
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-^RR - -^^R + v{R, r)v + fxiviV (2) 

where if) is the complex amplitude of the electric field 
in a nonlinear optics setting or the macroscopical wave- 
function of a Bose-Einstein condensate in an ultracold 
atomic gas setting. In our approach we focus on the latter 
framework, in which case V{R, T) is the external poten- 
tial. We will mainly consider the defocusing case a — 1. 
We note that Eqs. ([T]) and ([2| are not connected by a 
general transformation combining point-, gauge- or scale- 
transformations, which were discussed, e.g., in Refs. [13j . 
and M- 
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The variable R in Eq. ^ is strictly positive, so to 
simplify the treatment it is useful to apply the transfor- 
mation a; = In (i?) and t ~ which results in 



-V-x. + ^^(x,^)V' + IV'lV• 



(3) 



From now on we will consider the generalised equation 

1-04 = -ai{x,t)'4!^x + a2{x,t)ip + a3ix,t)\tljfip, (4) 

where ai_3 are some functions of x and t, and consider 
our case, Eq. as the following special case: 



o;2{x,t) = V{x,t), 
We now use the similarity ansatz [T5HT7] 



a^lx, t) = 1. 

(5) 



(6) 



where p G K"*" and {(p, 77} G K are some functions to be 
determined, and ^(77) is assumed to satisfy 



(7) 



where G is a constant and g and h are arbitrary functions 
of 77. By a simple transformation we can set (7(77) = 
without loss of generality. We have, for example, the 
following special cases: 



- 4>m + f70 + 20^ = 
with g = 0, /i = 77, G = 2 [18] and 

- 4>n-n + (j)^ = p,(j) 



(Painlcvc II) 



(Canonical NLS) 



with g = 0,h = — /i, G = 1. 

Substituting Eq. ^ in Eq. Q we obtain the following 
set of equations: 



(p2)t + 2ai(pVx)x = 0, 
r]t + 2(pxrixai = 0, 
—asp^ + Grjlai = 0, 

ai((/7^ - h-ql - — ) + q;2 + <(5t = 0. 
P 



(8a) 
(8b) 
(8c) 
(8d) 
(8e) 



Assuming Eq. ([5|, we obtain from Eqs. (8a)-(8d) the 
following solutions: 



e 3 



P = 
ip = 
77 = 



8ci{t) 
3 



2VGci(t) 



2x 

e 3 + C3, 



(9) 
(10) 

(11) 



where ci_2(i) are arbitrary functions of time and C3 is a 
constant. From Eq. (8e) we obtain an expression for the 
potential: 

1 , /.(77)e-^ 3 cl{t)+2c^{t)c^{t) 2, . ... 
"^=9^ + ^^- 16 W) 

(12) 

We will next consider the potentials and solutions that 
can be obtained with various choices of ci^2{t) and hirj). 



III. THE RING SOLITON-LIKE SOLUTIONS 

There is freedom in choosing the potential function 
Q!2. For example, the l/(9i?^) term can be removed if we 



choose h{rii) = and C3 = 0. Then from Eq. (12 1 we 



obtain a harmonic potential 

3 c2(t) + 2ci(i)ci(i) 



a2 



16 



clit) 



C2{t). (13) 



See Table |T] for a selection of choices for h{ri). 



A. The dark ring soliton-like solution, h{ri) = —p, 

This choice gives the canonical NLS (no external po- 
tential) so that from Eq. ([t]) we are now requiring (jjiv) 
to satisfy 



p,(f> = -0^^ -I- G(j) 



(14) 



where p > 0. 

Using Eq. (11), we can write down 0(77) = 
yZ/Ztanh [\/^('7 — 770)] , the kink solution of Eq. ( 14) with 
G = 1, as 



//i tanh 



/i3 / 63^ 63^0 

22 \^)~ ci(io) 



(15) 



re 3 ^° to match with 



where we have chosen C3 = 770 — 2c-i(to) ' 
the soliton centre xq at time t = to. 

Therefore, using Eqs. (|9|, ([To]) and wc have 
constructed the solution 



ci(T)i?3 



X tanh 



2^/2^' \ci(r) ci(To) 



(16) 



of the original radial GPE, Eq. ([2| (cr 1), with the 
potential 



V{R,T) = 



9i?2 cl{T)R3 

3 cj(r) + 2ci(r)ci(r) 

16 c?(T) 



(17) 



R^-C2iT). 
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TABLE I. A range of potentials available by suitable choices of h{ri). From Eq. ([7| we then have the equation for (t>{ri), which 
we must solve to construct the solutions '(l^{x,t). 

h{ri) C3 Potential Eq. to solve 

art -ID cji.jj 2G2c3(T) 



5 



2G2c3(T) V 4'? / 

1 -1_ -M 3 a;(T)+2c;i(T)ci(T) p2 ■ /^n , , , , ^ ,3 _ n 



Let us select ci(r) = 1 and C2(r) = -T. Then from 
Eq. ([l7| we get (see Fig. [l]) 



(18) 



and from Eq. ( 16 ) 



ip{R,T) = i?"3VAitanh 



3 



/^(i?3 -i?|) 



(19) 

where /i > (see Fig. [2]). 

We note that it is possible to rewrite the l/(9i?^) term 
in the potential as a fractional vorticity of ^ and include 
the full Laplacian. Then 



tp{R,T,9) = i?"sy^tanh 
solves 



2V2 



(20) 



iij^T = -^R.R - ^^B. - ^^ee + V{R)ij + IV-P^, (21) 



where 



V{R) = -M-T + 1- 
R3 



(22) 



B. The stability of the dark solution, h{rj) = — 

Direct propagation of 'ip{x,y,t) given by a two- 
dimensional extension of the GPE equation ([T]) in rect- 
angular coordinates y shows that the solution given 
in Eq. (19) is long-lived, until it decays into a vortex- 



antivortex necklace by the snake instability. The decay 
is much faster if instead of the potential in Eq. ( 18 ) we 
use e.g. 



ViR) 



A B 

R^^R' 



(23) 



DC 1 
> 



1 \ 




1 \ 
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FIG. 1. The potential (18 1 for a few values of ^. There 
is a minimum at finite R only for positive values of fj,. The 
potential approaches unity as 7? — > 00. 



where A and B are chosen so that V{R) and V{R) agree 
up to the first (linear) order when expanded at R = Rq. 

To get more insight on the stability we consider the 
Bogoliubov-de-Gennes equations for low-lying modes e [2 
131 IS] . In dimensionless natural units {jj, = h = 2m = 1) 
they are 




(24) 



where we have defined the Bogoliubov amplitudes u and 
V for which 



V'(r,T) 

by m 



^{'0o(r) + [u{r)e 



-ieT 



ieT 




Ug{R) 

MR), 



]}, (25) 



(26) 



representing a partial wave of angular momentum q rel- 
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FIG. 2. The ring dark soliton-like solution ( 19 1 with Ro 
and fi = 1. 
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ative to the condensate, and 



1 A 
RdR 



jL 

i?2 



1 



-2|V^oP, (27) 



and ipo{R, T) is given by Eq. (19 1 (with the substitution 
/i — !• A). If e > 0, w and v are orthogonal and normalised 
according to 



dr(|w(r)|2 - \v{r)\'^) = 1. 



(28) 



The instability time is given by Td — 1/Im(e), and if all 
the eigenvalues are real, the condensate is dynamically 
stable. The results of numerically solving Eq. (24) are 
shown in Fig. [sj There is only one (purely) imaginary 
eigenvalue, and the Bogoliubov amplitudes are localised 
at the radius Rs- Equation ( |24| ) is solved numerically 
using Lagrange functions, choosing a Laguerre discrete- 
variable representation [201 HI] , which is particularly effi- 
cient for radial geometry. The abscissas can be found us- 
ing Newton's method [22| , and the resulting linear eigen- 
value problem is solved by Hessenberg QR iteration [53] . 

1 occurs 



The minimum of the potential ( 18 ) with A 



at i?o = (1)^ w 0.44. Therefore, the soliton-like solu- 
tion ( jlQ] ) is stable if Rs is around (a healing length away 
from) the minimum of the potential. Another region of 
stability appears to be the hmit Rs ^ 00 (see Fig. |3]), 
although it is approached rather slowly. In this limit 
we obtain the one-dimensional Gross-Pitaevskii equation 
with a constant potential, which can be removed by re- 
defining the zero of energy. 

The case for Rs = 50.0, A = 1, and g = 8 is shown 
in more detail in Figs. |4]and[5] Because the ring suffers 
from the snake instability in general, and because there 
is only one mode with dynamical instability (per q) , they 
must be related. That this is indeed so is confirmed in 
Figs. |4] and [sj The amplitudes u and v corresponding to 
the imaginary eigenvalue are localised at the notch. Here 



the azimuthal symmetry was broken by choosing (261, 
that is, an overall phase of 0. 



FIG. 3. (Colour online.) Results of solving the Bogoliubov 
equations (24l for t he e xact solution (191 of Eq. ([2j| (a = 1) 



with the potential (18 1. Here A = 1. Shown are the decay 
times 1/Im(e), where e is the Bogoliubov eigenvalue. There 
is only a single imaginary eigenvalue per q and the corre- 
sponding amplitude is localised at the notch of the ring dark 
soliton-like solution. For low radii {Rs < 1.5003) the only 
unstable mode is g = 0. For radii 4.613 < Rs < 8.620 the 
primary decay channel is the quadrupole, or g = 2, mode, 
while all higher modes are dynamically stable. As the radius 
is increased the higher modes become unstable as well. The 
solution ( 19 1 is always formally unstable with respect to the 
q = mode, although the instability times are very long for 
large radii. 
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FIG. 4. (Colour online.) The amplitudes u and v correspond- 
ing to the only imaginary Bogoliubov mode, and the wave- 
function ipo for Rs = 50.0, A = 1, g = 8, and e = at T = 0. 
The inset shows the Bogoliubov spectrum e. The amplitudes 
It and V are localised around the notch at Rs = 50.0. We 
have chosen one possible normalisation satisfying (|28[|. 



C. The bright ring soliton-like solution, h{ri) 



When G — —1, Eq. (14 1 has also a bright soliton- 
like solution (j){r]) = \/2-^ — /isech [-y 



m)] (m< 0). 

-1 so nothing 



Similarly as in Sec. Ill A (now 
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FIG. 6. The ring bright soliton-like solution (301 with 7?o = 15 
and ^ = —1. 



with the potential (18), where fi < (see Fig. [6|), of the 
original radial GPE, Eq. ([2]) (cr = -1). 



IV. PAINLEVE II, h{r]) = r? 
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FIG. 5. (Colour online.) The instability mode is directly re- 
lated to the snake instability, having typical behaviour around 
the ring (for q = 8). |(a)| The amplitudes u and v of Fig. [4] 
added to the wavefunction tpo at 9 — 0. The disturbance is 
towards increasing radius. The ins et sh ows a magnification 
of the area inside the red rectangle, (b) At 6 — it/8, the dis- 
turbance is now towards decreasing radius, (c) At S = 27r/8, 
the disturbance is again towards increasing radius. |(d)| At 
9 = 37r/16, between the previous extremes, the disturbance 
does nothing, (e) At 6 = 77r/32, which is not an extremal 



angle, the disturbance is also between extremes. 



changes), we arrive at the solutions 



ci(r)i?3 



X sech 



ci{T) c,{To) 



(29) 



with the potential (17) (/i < 0) and 

"3 



ipiR,T) = i?"5y^sech 



-^liR--~R§) 



(30) 



We are now requiring 0(77) to satisfy 
- 0«« + ry<^ + 24>^ = 0. 



(31) 



This means 3(77) = 0, h{rj) = rj, and G = 2. As we can 
see m Table |l) Eq. ^ arises as the similarity reduction 
of Eq. ([2]) (cr = 1) with 



V{R,T) 



1 



9i?2 16 
3 



3 c?(7^)+2ci(T)ci(r)^2 



c?(T) 



(32) 



2ic?(T) 



-C2iT). 



This is similar to what happens when the modified KdV 
equation 



Vt -6v + Vxxx = 

is transformed with a similarity reduction 

, , w(z) X 

v(x,t) = 7-, z = r 

(3^)5' (3<)5 

such that w(z) satisfies Pn: 

Wzz = + zw + a. 



(33) 



(34) 



(35) 



where a is a constant of integration [21] . 

Interestingly, we can eliminate time dependence from 
the potential ( 32 ) without sacrificing nontrivial time de- 



pendence in the solution, unlike what happened in Sec- 
tions III A and |III C| by choosing 



ci(T) = C4Cosh [2Vfc(r-f Cg)]^ 



(36) 
(37) 
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where k ^ Q, C4, C5, and cg are constants. 
Therefore, 



V^(i?,r)-p(i?,T)e»'^(«'^V(r7(i?,r)), 



where 

p{R,T) 



C4Cosh \2Vk{T + C5)]i 



■tanh [2VA?(T + C5)]i?2 + C7 



c&{T + C5) 



3 tanh [2^/fc(^ + C5)] 



^/fc23c3 



2ic4COsh [2yfc(r + C5)]; 
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(38) 
(39) 

(40) 
(41) 



and where c-j is a constant, solves Eq. Q (cr = 1) with 
the potential ( 32 1 



Vm = ^,-kR 



C6, 



if 0(77) solves Eq. ( |3l[ ), the second Painleve equation. 
If A: < 0, we obtain from Eq. (36) 



ci(r) = C4COS [2v/|fc|(T + C5) 



Then Eq. (381 with 



C4COS [2^|fc|(T + C5) 



2 5 
3 



tan[2v/|fc|(T + C5)]i?'^ + C7 



■C6(T + C5) + 



p(i?,T) 
^(i?,T) 

r;(i?,T) = ^= 

23c4Cos[2^(T + C5)]3 

solves Eq. ([2]) with the binding trap potential 



3tan[2v/|fc|(T + C5)] 
k\2icl 

-Ri 



V{R) = 



9i?2 



|fc|i?2 



(42) 
a. 

(43) 
(44) 

(45) 
(46) 

(47) 



Any nontrivial real solution of Eq. ( 31 1 with the bound 



ary condition (j>{r]) — >■ as 77 — ^ 00 is asymptotic to vK\{r]) 
for some nonzero real constant v, where Ai denotes the 
Airy function |5S]. The choice = 1 corresponds to the 
Hastings-McLeod solution [26], and the choice ji'l < 1 



to the Segur-Ablowitz solution [27]. Since by Eqs. (41 1 
and (46 ) we have > 0, the solution is well approximated 



by taking the Airy function for (j){rf): 
(j){r]) = i'Ai{r]). 



(48) 



From Eqs. (|38|, (|44|, (|46|, and (|48]) it follows that 
R\ij\'^ for aU R as 



T 



C5, 



(49) 



where I € Z. We have made density plots corresponding 
to v = 0.9 and C4 = 3/25 for various values of k (see 
Fig. [7]) . Note that unlike for the dark soliton-like solu- 
tions, here the solution is actually dynamical and evolv- 
ing in time, describing a scattering (one-time or periodic) 
of the bright solution from the central potential. 

The limiting potential with fc = can be obtained by 
choosing 



Cl(T) =C4(r + C5)s, 
C2{T) ^ 



25ci(T)3 



C6. 



Now Eq. (|38]) with 

R-i 

C4(T + C5) 



cp{R,T) 



1 



T + C5 
C7, 



i?2 



3V2 
8^ 



R 

T + c, 



solves Eq. ([2| ((7 = 1) with the potential 



(50) 
(51) 



(52) 
(53) 

(54) 
(55) 



See Fig.[7][(f) 



We note that the k = case is also obtained by choos- 
ing ci(T) = const, and C2(T) — cq, which leads 

22 cl 

to 



(56) 
(57) 

(58) 
(59) 




cejT + cr, 



In this case the equation may be called a "nonlinear 
Bessel" equation. 



V. SUMMARY 

We have obtained exact localised solutions of the radial 
Gross-Pitaevskii equation describing cylindrically sym- 
metric systems. We have concentrated on solutions that 
may have implications for studies of ring solitons in quan- 
tum gases, but it should be noted that still other solu- 
tions are possible with different choices of h{ri) and the 
integration "constants" ci_2{T). 
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FIG. 7. (Colour online.) Density plots of -R|-0p as given by Eq. (38 1. In all of the plots u — 0.9 and C4 = 3/23. In (a) and |(b)| 
the radius of the centre of mass is oscillating and the condensate is breathing between a point and an extended state. The 
magnitude of k determines the period of the oscillations. In (c) (d) and (e) there occurs a reflection at the repulsive core 



potential around i? = of an incoming density from infinity. The magnitude of k determines the spread and velocity of the 
incoming and outgoing rings and the duration of the reflection dynamics. C5 determines when the reflection occurs in time. 
In |(f)| is shown the limit fc — )■ of potentials ( |42[ ) and ( |47[ l, the limiting behaviour of oscillations as the period goes to inflnity. 
Note that the colouring is not to scale. 



The form of the potentials corresponding to exact so- 
lutions suggests that confinement of cold atoms in traps 
which are repulsive at the core should be studied further, 
especially in the context of ring soliton stability. Previ- 
ous work and our own numerical simulations have shown 
that in general the ring dark soliton is ultimately de- 
stroyed by the snaking instability, and the ring collapses 
into a vortex necklace . Having a solution with a po- 
tential that has both a repulsive core, and a minimum at 
finite R, however, suggests that e.g. toroidal traps (also 
known as ring traps) might provide better stability for 
the ring dark soliton, and results of numerical investiga- 
tions testing this idea will be reported in a later paper. 
Such potentials are experimentally feasible [28 30 . We 
will also discuss elsewhere the possibility of creating dark 
solitons in toroidal traps using time-averaged potentials 
and/or adiabatic passage techniques [5m34j . 

An additional aspect to such studies comes from the 
observation that some of our solutions allow one to choose 
freely the location of the dark or bright solution (Rq), as 



this location does not depend on the parameters of the 
potential. Thus the long-time survival of the solution is 
not connected to locating the soliton-like structure at the 
minimum of some confining potential. It means that one 
may study ring soliton dynamics before the decay into 
vortices takes place, and look for similarities to such be- 
havior in one-dimensional systems [35]. Another feature 
is to consider the case of a bright ring soliton in such 
traps, as it has not been studied very much in the past. 
Such considerations form the starting point for further 
work on ring solitons and soliton-like structures, includ- 
ing their stability and dynamics. 
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